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Abstract 

We study the holographic dual of a (2 -|- 1)-dimensional s-wave superfluid that 
breaks an abelian U{1) x U{1) global symmetry group to the diagonal U{l)y. 
The model is inspired by Sen’s tachyonic action, and the operator that condenses 
transforms in the bifundamental representation of the symmetry group. We focus 
on two configurations: the first one describes a marginal operator, and the phase 
diagram at finite temperature contains a first or a second order phase transition, 
depending on the parameters that determine the theory. In the second model 
the operator is relevant and the finite temperature transitions are always second 
order. In the latter case the conductivity for the current associated to the broken 
symmetry shows quasiparticle excitations at low temperatures, with mass given 
by the width of the superconducting gap. The suppression of spectral weight at 
low frequencies is also observed in the conductivity associated to the conserved 
symmetry, for which the DC value decreases as the temperature is reduced. 
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1 Introduction 

The realization of solutions with scalar hair in asymptotically AdS (aAdS) black holes in [1] 
prompted the study of strongly coupled superfluid systems via the holographic gauge/gravity 
duality. The non trivial profile of the scalar field corresponds in the field theory to the 
condensation of an operator. If, in addition, the scalar is charged the model contains the 
breaking of the global symmetry dual to the gauge symmetry in the bulk, hence superfluidity. 

The seminal work [2] considered an Einstein-Maxwell-Higgs (EMH) action, with the scalar 
field coupled minimally to the U{1) Maxwell field (see [3, 4] for comprehensive reviews). 
Further models with order parameters of different spin (s- [5], p- [6, 7], and d-wave [8]) have 
been developed in the last years. And while the initial setups in refs. [3, 4] are bottom-up 
models with the minimal field content to describe the superfluid phase transition, microscopic 
embeddings have been proposed in the framework of type IIB string theory [9], M-theory [10], 
and D-brane models [7]. 

In this work we draw intuition from string theory to modify the minimal setup of [2] and 
construct a model which realizes superfluidity as the spontaneous breaking of the 17(1) x 17(1) 
flavor symmetry supported by a D3-D3 brane pair in an asymptotically AdS^ black hole 
geometry {aAdS^ BH)^. Effective brane-antibrane actions derived from string theory [11] were 
first used in [12] to realize chiral symmetry breaking through tachyon condensation in a stack 
of overlapping branes and antibranes (see also [13] for further developments in this direction). 
More precisely, those effective actions contain, in addition to the worldvolume gauge fields, 
a complex scalar, the tachyon, which is the lightest open string mode extending between the 
brane and the antibrane. In the models [12]-[14] the infrared region of a confining geometry 
makes the tachyon condense and consequently break U{Nf)L x U{Nf)ji —)• U{Nf)y. 

^See [7] for models of holographic superfiuids based on two overlapping probe branes. 
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We consider a setup corresponding to a single spacetime filling D3-D3 pair in an aAdS 4 
BH, and switch on a chemical potential along the axial U{1)a C [/(I) x U{1) living on the 
worldvolume of the branes. We couple the branes to gravity and take into account the 
backreaction, studying the system as the temperature is varied. While at large temperature 
the setup is described by a solution with vanishing tachyon, below a certain temperature 
the dominant solution contains a non trivial tachyon that realizes the spontaneous breaking 
U{1) X C/(l) — U{l)v, needed for describing a superfluid phase transition. 

This model possesses some characteristics that add novelty and interest to the analysis 
carried out in this article. First, unlike most of the examples in the literature, in this setup 
the dynamics of both the scalar and the gauge field are governed by a non-linear DBI term.^ 
Second, our setup is dual to a state with a non-zero density of matter transforming in the 
fundamental representation of the gauge group, and the tachyon (corresponding to the order 
parameter of the superfluid phase transition) is dual to a bifundamental operator. Third, our 
model allows not only to tune the charge and the dimension of the operator that condenses 
(dual to the charge and mass of the tachyon), but also to vary the strength of the backreaction 
of the matter sector of the action. As a result a rich phase diagram is found. Finally, the 
presence of a second U{1), the U{l)v not coupled to the tachyon, opens up the possibility of 
studying unbalanced superconductors; as in [17]-[19] one could switch on a chemical potential 
along the unbroken U{l)v and interpret the system as an unbalanced mixture. 

The paper is organized as follows. In section 2 we present the action of our model and de¬ 
rive the equations of motion. Section 3 describes an analytic solution with vanishing tachyon, 
dual to the normal phase of the system. There it is also shown that for some values of the 
parameters in the action, corresponding to different field theories on the boundary, the system 
is unstable towards condensation of the tachyon. To end that section we discuss how to set up 
the numeric integration of the condensed phase. In section 4 we analyze the condensed phase 
for the case in which the scalar is dual to a marginal operator. We are able to identify the 
IR geometry of the zero temperature solution (AdS^) and to construct the phase diagram of 
the system. Both first and second order phase transitions are shown to occur, depending on 
the values of the parameters in the action. The case of a scalar dual to a relevant operator of 
dimension A = 2 is discussed in section 5. In this case the phase diagram presents only second 
order phase transitions. Section 6 is dedicated to the study of the conductivity. Numerical 
results for the two theories studied in previous sections are presented there. We conclude in 
section 7 with a discussion of our results and possible continuations of this work. 

2 Setup 

In this section we present the model under consideration; Einstein gravity coupled to the 
tachyonic DBI action describing a spacetime filling D3-D3 brane pair. The starting point is 
then Sen’s tachyonic DBI action [11] (see also [20]-[22]) 

5 =-y'd^xR(|r|) (v^^d^07^+\/^d^tA^) , (2.1) 

with V a potential to be determined below and matrices 

-k , (2.2) 

^In [15, 16] holographic superfluids where the dynamics of the gauge field is described by a Born-Infeld 
term were considered for the first time. 
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where is the metric, {A^) the U{1) gauge field in the worldvolume of the (anti-)D- 
brane with field strength = dA^, and the complex scalar r is the tachyon.^ The tachyon 
is charged under the axial U{1)a subgroup of the original 17(1) x 17(1), and neutral under 
U{l)v, 

D^ = d^ + iq{Al-Al) . (2.3) 

We can now make explicit use of the combinations 

A^ = A^ + A\ A^ = A^-A^, (2.4) 

defining the vectorial 17(l)y, and axial 17(1 )a gauge fields respectively. In terms of these one 
can write 


— G'/xi/ + + F^^ + , (2.5a) 

= G^u + , (2.5b) 

D^ = d^ + iqA^ . (2.5c) 

In the equations of motions obtained from (2.1) it is consistent to set A^ = 0, hence the 
action (2.1) can be truncated to 

S = -2j d^xV{\T\)^-det[G^, + F^, + D^^TD;^ , ( 2 . 6 ) 


where the only gauge field is now the axial combination of the original 17(1) x 17(1) symmetry. 

For our setup to be dual to a 3-dimensional system at finite temperature and chemical 
potential, the D3-D3 brane pair described by (2.6) will be embedded in an asymptotically 
AdS^ black hole geometry. At low values of the temperature (measured with respect to the 
axial chemical potential) the expectation is that the backreaction of the charged branes on 
the geometry becomes important [23]. To account for this fact we couple (2.6) to Einstein 
gravity, which fixes the action of our model to be 


S 






V^(R-2A 


-TbV{\T\)^J-det[P[G]^u + F^l,+ , 


(2.7) 


In this action 2A = —djl? is a cosmological constant and F[G] stands for the pullback of the 
metric onto the branes worldvolume (from now on we consider the static gauge = G^y). 

Notice that we have written explicitly the tension of the brane in the Lagrangian. In our 
bottom-up model we will consider it a tunable parameter which regulates the amount of 
backreaction by the branes. As for the tachyon potential we will take it to be 


I7(r) = exp (m^ jrj^) , (2.8) 

which for small values of r behaves as in [21], namely 

l + mVV + 0(r*r)2 , (2.9) 


where m is then the mass of the tachyon, which is related to the dimension of the operator 
dual to r in the standard way in holography, A(A — 3) = wiAl?^ with L the radius of the 
asymptotic AdS^ spacetime. 

®In the notation of [21] we have set ttq' = 1 and the transverse scalars to zero. We follow conventions for 
(anti)symmetrized indices with no factors of 2, i.e., M^ab) = Mab + Mi,a, etc. 
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Notice that in this bottom-up approach we will take to be an order one parameter. 

As a consequence this model cannot be constructed from string theory without including a' 
corrections [13].^ 

Finally, notice that the radii appearing in and in the cosmological constant, 2A = 

—6/L^, are not the same. To relate them consider = r = 0 in (2.7), then the cosmological 
constant is shifted 


- 2A ^ -2A 


n = -2A = ;^ 


( 2 . 10 ) 


In the solutions of (2.7) it is the shifted cosmological constant A the one that controls the 
AdS^ radius near the boundary. In the following we will work in terms of the AdS^ radius L 
defined as in (2.10). 


2.1 Ansatz and equations of motion 

We will be interested in static solutions with rotational and translational invariance, conse¬ 
quently 

G^udx^dx’^ = Gtt{r)dt^ + Grr{r)dr‘^ Gxx{f’)dx^ , ( 2 - 11 ) 

with negative Gtt- Moreover, we will set the phase of the tachyon to zero, and consider 
solutions where the only non-vanishing component of the gauge field is the temporal one: 

r* = r = r(r), At = At{r). (2.12) 


Let us now write the equations of motion for this ansatz. First, it is useful to define 


W = {G + F + = M-^ , 


with components written with two risen indices: In particular we need 


(2.13) 


n = 

w** = 

yyhi] = 


Vdet Ad = 


[Af+ (G« + 2g%V)(G.. + 2T'2)] , 
Grr + 2 t '^ _ 2q^A^ t"^ + Gtt 

n ’ “ n 

— 2 / 1 ' 

= G-i , 

Gxx n 

^Gtt + 2q^AfT^VGrr + 2r'^ ■ 


(2.14a) 

(2.14b) 

(2.14c) 

(2.14d) 


Then, the equations of motion read 


Rfiu — -|- ( A — 

1 


) = 0 . 


V{t)\/— det A4 
1 

V det A4 


dr 

dr 


F(r)V-detA4>v[^*'l = 0 , 


F(r)V-detA4 W'''' r' 


drVjr) 

2V{t) 


+ W^^q^Ajr ) = 0 . 


(2.15a) 

(2.15b) 

(2.15c) 


^Moreover, according to the original conjecture for the tachyonic action, the potential evaluated at its 
minimum should vanish [20]. Although this occurs for the confining geometries considered in [12]-[14], it will 
not be true for our superfluid solutions. More differences with a top-down approach consist of the absence of 
a dilaton or Wess-Zumino terms. 
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Next,we partially fix diffeomorphism invariance by choosing as metric ansatz 


ds'^ = -g{r)e + ^dx^ + , 

g[r) 


(2.16) 


and the equations of motion (2.15) then reduce to the following set of two second order and 
two first order differential equations 


, a / Tb\ nrV(T) 
+ - + ( -3-y ) + 

n 

^ g^S 


2S 

q^n 


(1 + 25 /^) = 0 , 


^ r l/(r) H- re^V (r) A'f = 0 , 


// 9 , , /N /3 , / 2 , 1 

(rgx - 4Q-rg)T'^+ I --^-h 


r — 


rQ 


rQ 


2 \g ' QJ 2Q ' 


^ e'^ Af V' (t) ^ q‘^ Aj T drV{T) ^ f 2q‘^ Aj t drV{T) 


2QV{t) 


gQ 


2gV{T) 


Q 


Vir) 


t'2 = 0 , 


A': + 


2(5 - Q)'r' 2 o' /I '^\,gX'9r' g , , 


Qt r 

2g2 At 2{g - Q) Af g{Q - g) 


+ 


+ 


Af -4q^ AtT^r''^ = Q , 


2 _/2 


q^ rQA^r'^ ' 


g AtQ 

where we have defined 

Q = g-2q^ Aj , 

S = \j\ — eA Af + 2 5 r'2 — [2q^ / g) AlT‘^{l + 2g r'^) . 


(2.17a) 

(2.17b) 

(2.17c) 

(2.17d) 


(2.18a) 

(2.18b) 


As we will now see, in some limits these equations reduce to those of well know systems 
as the holographic superconductors of [5]. 


Linear limit 

It is possible to recover the original holographic superconductor setup of [5] in a limit in which 
the DBI term in the action is linearized. This linearization is achieved by taking Tf, —>• 00, 
with 

^ ' ~ fixed , -\/7l A —)■ A fixed , \/^b t —)■ r fixed . 




q 


(2.19) 


In this limit the action (2.7) reduces to 


5 = 




- D^fD>^f - + 0{T^^) 


( 2 . 20 ) 


where we have used the relation (2.10). Indeed this action is that of the minimal model for 
holographic superconductors constructed in [5]. 


Probe brane limit 

In the opposite limit, Tt, —)• 0, the equations of motion for the metric decouple from the equa¬ 
tions of motion for the 1/(1) gauge field and the tachyon. Therefore, in this limit our system 
can be interpreted as a spacetime filling D3-D3 pair in the Schwarzschild-Ad £'4 spacetime, 
with the action describing the dynamics of the branes given by the DBI part of (2.7) alone. 
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Unbalanced superconductors 

In [18], a model realizing unbalanced holographic superconductors was constructed. In unbal¬ 
anced superconductors two different fermionic species with unbalanced populations contribute 
to the formation of superconducting states. This imbalance, or chemical potential mismatch, 
can be described holographically by means of a second 17(1) under which the order parameter 
is neutral [18] (see also [17]). 

In our model, in addition to the chemical potential along the U (1)^ directly coupled to the 
tachyon (see eqs. (2. 4), (2. 5c)), one could switch on a chemical potential along the remaining 
U{l)v- The tachyon is not charged under this second 17(1), and thus following [18] one could 
treat it as dual to a chemical potential imbalance. In [18] the two 17(l)s did not interact 
directly with each other, but only through their backreaction on the geometry. This will be 
different in the present system due to the non-linearity of the DBI action. Hence, it would be 
interesting to study the phenomenology of this model when switching on a chemical potential 
along U{l)y. We do not pursue this direction in this paper, but leave it instead for future 
works. 


3 Solutions 

In this section we present an analytic solution of the equations of motion which is dual to 
the normal phase of the system (for which the charged scalar vanishes). Next, we show that 
at low temperature this solution is unstable towards condensation of the scalar. We study 
the dependence of this instability on the parameters of the model, which will be crucial for 
the characterization of the phase diagram in the following sections. Finally, we study the 
boundary conditions which will allow us to numerically construct the solutions with a non 
trivial scalar [i.e. those dual to the condensed phase). 


3.1 Normal phase solution 

A charged black hole solution to (2.15) with r = 0, corresponding to the normal phase of the 
dual theory, was found in [24] 


g{r) 


L2 


('•i? 


nL'^\p\- 


2F1 


( 11.5. 1 

2’ 4’ 4’ ) 

2^2 



(3.1a) 


X = 0 , 

P j—Y A -J 


(3.1b) 

(3.1c) 


where p is a constant of integration related to the charge density of the black hole and kh 
is a constant of integration that can be translated into the radius of the horizon, r_|_, by the 
condition p(r_|_) = 0 (with r_|_ the largest root). 

The thermodynamics of this solution was analyzed in [24, 25], with temperature, entropy 
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density, chemical potential, charge density and pressure density given respectively by 

1,2 


T = 


r+ 

dvrL^ 


1 - 


3 + TbL^ 




8tt 


271 ri 

’ 

/i = p - 2F1 

r+ 

n 


1 15 
2’4’4’ 


r4 
2 ^ 
-P — 
rl 


P = 


1 1 


-k,,-^\pfAr{-3/A)r{5/A) 


(3.2a) 

(3.2b) 

(3.2c) 

(3.2d) 

(3.2e) 


3.2 Instabilities towards scalar condensation 

We now want to study possible instabilities of the normal phase solution (3.1) towards the 
condensation of the tachyon. Let us then consider the fluctuation of the scalar field r around 
that solution. This perturbation is described by the DBI part of the action (2.7) with the 
metric and U{1) field as in (3.1). It was shown in [26] that to describe fields appearing in the 
DBI action it is convenient to work in terms of the open string metric, defined as the 
inverse of the symmetric part of the LV(r = 0) matrix (2.13) (see [20] for similar comments 
in the context of Sen’s effective tachyonic action). Equivalently, 

" =G^,-{F- G-i 




F) 


flU 


(3.3) 


A second quantity that is important for the description of the perturbation is the effective 
running coupling 

= 7-det(G + F) ’ 

that compensates the change of volume density from the original description to the one in 
terms of the open string metric. For the solution (3.1) 


13 = 




Dfj_udx^dx’^ = ( —gdt^ + + ^dx^ , (3.5) 


with g{r) given in (3.1a). The equation of motion for the perturbation of the scalar field then 
reads 




■■di, 


'—s 




- [rrA + q^Afs^^) r = 0 . 


(3.6) 


Close to the boundary At 
given by m?. 


p, and (3.6) describes a free scalar in AdS 4 with mass squared 


Instability at extremality from the AdS 2 throat 

At zero temperature the horizon is a double root of the blackening function (3.1a). Therefore, 
close to the extremal horizon (at r = tq) both /3 and approach a constant. Hence, in 

the zero temperature limit the equation of motion (3.6) describes a scalar with effective mass 


rrigg = lim (m^ + q^Afs^^) = rrA — 2q 


r—>-ro 


, 6 + Tfe 


(3.7) 
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Figure 1: Critical value of the charge as given by = “1/4 (solid line) and by the 

onset of the instability in the aAdS^ region (dashed line) for two values of the scalar mass. 
The shadowed area is unstable towards condensation of the scalar field. 


in the near-horizon geometry AdS 2 x with AdS 2 radius 


r2 

^AdS2 




2 


6 (6 + rbL2)(6 + 2rfeL2) 


(3.8) 


As was first shown in [27, 28], there is an instability when the effective mass violates the 
Breitenlohner-Freedman (BF) bound of AdS 2 - ^ “1/4- The violation of the BF 

bound defines a critical charge (/crit; dependent on the mass of the scalar field and the tension 
parameter T^, above which the scalar r is unstable. In figure 1 we plot this critical value 
as a function of the tension Tb both for the theory with m? = —2 and m? = 0. This 
results in a curve (solid line) above which the normal phase described by eq. (3.1) is unstable 
towards condensation of the scalar. 

From the arguments in [29, 30] the instability caused by the violation of the BF bound of 
the AdS 2 near-horizon region gives rise to a continuous phase transition. In particular, near 
the extremal horizon tq the equation of motion (3.6) for the scalar field takes the form 


dr {{r - ro)V') - mlsL\^s 2 ^ = 0 , 
which, for < —1/4 has as solution 


r = 


y/r - ro 


sm 





(3.9) 


(3.10) 


where ruy is a UV scale at which the AdS 2 region has to be corrected to recover the AdS/^ 
asymptotics of the full solution, and Cr is a normalization. 

As explained in [29, 30], equation (3.10) implies that there is a quantum phase transition 
if we allow the tension of the brane to vary from a BF bound-violating to a non-BF bound- 
violating value. For tensions slightly into the unstable region (he., < —1/4) the 

critical temperature at which the phase transition takes place and the value of the condensate 


















at vanishing temperature are given by 


Tc ~ Auu exp 





— TT 

, {Or) ~ Air exp 

— TT 

1 

+ 

_1 



(3.11) 


with Aiij (Auv) an IR (UV) scale. 


Instability at extremality from the boundary of AdS^ 

Let us now look again at figure 1, which illustrates the conclusions of the stability analysis 
carried out in the previous subsection. That analysis established that for points in the {Ti,,q) 
plane above the critical curve (solid line), the normal phase solution (3.1) 

is unstable at small temperatures. For points below that curve the BF bound is not violated. 
In the linear model of holographic superfluidity [27] this region is stable and dominates the 
large temperature part of the phase diagram. However, in the present case a careful analysis 
shows that this no longer holds; for low values of the tension there is a region below the 
critical curve in which the AdS 2 BF bound is not violated, and yet there is an instability of 
the scalar, indicating the preference for a condensed phase. This is the region between the 
dashed and the solid lines in figure 1. 

To understand this new instability let us focus first on the probe approximation, Th = 0, 
at zero temperature. The system becomes that of a probe brane in AdS 4 ,, with an electric 
field on the worldvolume of the branes turned on. The open string metric describing the 
fluctuation of the scalar governs the form of the equation of motion for r. Near the origin 
this equation becomes 

,! 2 , L? 

T"+-T' + -^r = 0 , (3.12) 

which is the equation of a scalar in AdS 2 with effective mass Lp' = —q^ Lp (in units of 
the AdS 2 radius).^ The solution of this equation reads 


Ct 

r = —p sm 

Vr 



log 


r 

rvv 


(3.13) 


As one can see, there is a minimum value of the charge for this solution to exist, be real, 
and well behaved at the origin. Consequently, at —)• 0 there is an instability of the normal 
phase solution for 

(3.14) 

independently of the mass of the scalar (and the dimensionality of the theory). Interestingly, 
as we explain below, we have found that this extremal instability for qL > 1/2 extends 
to small values of Tf,. However, at non-zero T^, the instability is not anymore caused by a 
violation of the BF bound in the AdS 2 throat. It must then be related to the boundary 
conditions associated to the UV of the throat, i.e., the region where the throat opens up to 
the asymptotically Ad54 geometry. In particular this is an effect of the non-linearities of our 
model, since it is not present in setups like [27]. 

®One should keep in mind that, although effectively the scalar r behaves near the horizon as if it were in 
an AdS 2 geometry, the system corresponds to a probe in AdSi. 
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These generic boundary conditions allow for the unstable modes always present in AdS 2 to 
survive in the spectrum. A similar scenario is given by the instabilities studied in a different 
setup in [31], where the authors allow mixed boundary conditions for the scalar field, as 
opposed to the Dirichlet boundary conditions that we impose in the aAdS^ boundary in the 
present paper. 

We have been able to find the values of the charge at which, for a given tension T^, 
the T = 0 instability we have just described occurs. The result is given by the dashed 
line in hgure 1. This line is constructed by hnding the normalizable solution to (3.6) at 
several small temperatures and extrapolating down to T = 0. As we will see later, at very 
small temperatures the phase transitions are always second order, which guarantees that this 
prescription to find the dashed line in figure 1 gives the correct result.® 


3.3 IR boundary conditions of the condensed phase 

Once we have shown that the normal phase solution is unstable towards condensation of 
the tachyon the next step is to find the solutions with non trivial scalar, which are dual to 
the condensed phase of the system. To construct such a solution we have to integrate the 
equations of motion numerically. We can make use of the following scaling symmetries: 


t r \ g 

At 

L 

g/3 g7 1 g27 g2(a+7) 


g/3+7 


We use the y-scaling to impose that the metric approaches AdS^ with unit radius when 
r ^ oo 7 There are two remaining scaling symmetries which leave L invariant 


t ^ 6°^ t , X ^ X + 2q: , At ^ e ^ At , (3.15a) 

X* ^ x* , r ^ e~^r , g g , x ^ X “ 2/3 . (3.15b) 


We use (3.15b) to fix the position of the horizon at r_|_ = 1. The symmetry (3.15a) will be 
used once the numeric integration has been performed to relate the periodicity of time in the 
UV with the temperature. 

The asymptotic behavior of the fields in the IR is given by a series expansion whose leading 
terms at finite temperature are 


gr = e^+/^47rT(r — r+) + • • • , (3.16a) 

At = o+(r - r+) H-, (3.16b) 

X = X+ + --- , (3.16c) 

T = T+ + e -^+7 ^ ^ ^ . (3.16^) 

AttT r+' 


where x+) and t+ are unspecified constants and 


T = 


g'e 2 


dvr 


= e 


^1 - a^ex+ (6 + n) - 


r=r+ 


2./I - aiex+ 


(3.17) 


®If the transitions were first order we would have found the value at which metastable solutions cease to 
exist. 

^ We do so in the remaining of this paper, although we restore the AdSn radius for the dimensionless 
combinations Tt , q L and rr7 for presentation purposes. 
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is related to the temperature of the configuration. The reason why T is not readily identified 
with the temperature is that the appropriate boundary condition is given by the periodicity 
of the euclidean time on the boundary. Actually, near the UV we impose that g ^ and 
X —>• 0 to calculate the temperature, but this cannot be done until the full numerical solution 
is obtained. 

The strategy we follow to calculate the numerical solutions is as follows. From eq. (3.16) 
we have three free parameters in the IR: x+j aiid r+. While on the boundary we want to 
impose two boundary conditions: the vanishing of the non-normalizable mode of the scalar, 
and X 0. We will fix x+ = 1 from the beginning, and make use of the scaling symmetry 
(3.15a) once the numeric solution is obtained to ensure that x 0 at the boundary. Then, 
we expect a family of solutions given by a relation between t+ and a+ that ensures the correct 
boundary condition for r at infinity. The resulting one-parameter family of solutions will give 
us the value of the condensate as a function of the temperature of the system. 


3.4 Free energy and VEVs 

Before presenting our solutions let us explain how we can detect a phase transition and 
calculate the condensate of the scalar operator dual to r. 

When two or more solutions coexist at the same value of temperature and chemical po¬ 
tential we should determine which one is thermodynamically favored by comparing their free 
energies. This quantity is obtained holographically from the evaluation of the euclidean on- 
shell action on the boundary. From equation (2.15a) it follows that the euclidean on-shell 
action is 


^on-shell — 


/% R 2 


dry/gR^^ 







rA 



(3.18) 


where we have integrated up to a cutoff r\, V 2 is the spatial volume of and = 1/T 
the period of the euclidean time in the UV. The lower limit of this expression vanishes since 
g ^ 0 there, whereas the upper limit diverges when r\ ^ 00 , and must be regulated with a 
set of counterterms evaluated at the surface r = r\ 


/3r^2 


^ counterterms — 




K - 2 - — 
2 


r=rA 


(3.19) 


where 7 is the induced metric on the constant-radius slice, and K = Kab'y°'^ with Kab = (a^b) 

the extrinsic curvature and Ub the components of an outward-pointing unit vector orthogonal 
to the constant-radius slice 

Kah = dia.gi^e~^y/g{g'- g^) ,ry/g,ry/i^ , K ^ ^ 


The hrst term in (3.19) corresponds to the Gibbons-Hawking term, and is needed for the 
definition of the variational problem in gravity, the second term is a volume counterterm that 
cancels the divergent behavior of the asymptotic Ad 54 spacetime [32], and the last term is 
needed to cancel a divergence given by the scalar field for some values of its mass.® A term 
proportional to would produce a finite contribution to the free energy, but since the 

original action is even in the scalar we do not include it. 

®And in particular when it is the leading term near the boundary the one that is kept fixed in the variational 
problem. 
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The internal energy and pressure can be obtained as conserved charges from the Brown- 
York tensor [32], and in our case reduce to the tt and xx components of 

= + . (3.21) 

Finally, from the action (2.7) complemented with the counterterm (3.19) we can calculate 
the variation of the on-shell action with respect to the non-normalizable mode of the scalar 
field, tq. This variation is dual to the vacuum expectation value (VEV) of the corresponding 
operator 


(Or) = T 


dinn-: 


on-shell 

6to 


= -T 


F2 (E(r) Vd^t^ Wr' + 2V7 


6t 

Sto 


(3.22) 


In the following sections we will write down explicit versions of this expression and that of 
the free energy for the two cases analyzed in this paper (namely the m?L?‘ = —2 and w? = 0 
theories). 


4 Phase diagram for the model with a marginal operator 


In this section we discuss the case in which the scalar field is dual to a marginal operator 
with conformal dimension A = 3. Holographically this is obtained by taking the mass of the 
scalar to zero, rm? = 0. Then, the potential given in eq. (2.8) is just V = 1. 

In order to read the physical observables of the dual theory we need to compute the 
boundary expansions of the different helds in our setup. For a massless scalar the two leading 
terms of the independent solutions to the scalar equation near the boundary (r —)• oo) are 

r~roH --• (4-1) 

The source of the dual scalar operator is determined by tq and the VEV is given in terms of 
T 3 . We are interested in the normalizable solution with zero source, tq = 0, and we impose 
this in the following. 

With To = 0 the remaining equations of motion can be solved for large r as 


g r 


' £ 

1 H- n + 


^3Tbr| 

X^Xo + —^ + --- , 

A Q 

~ /r H-^-. 

r 


(4.2a) 

(4.2b) 

(4.2c) 


The free energy (H = -T (lon-sheii +-^counterterm)), internal energy (E), and pressure (P) 
obtained from these UV (r —)• 00 ) expansions, and the relations (3.18) - (3.20) are given by 

p-Xo/2 T? 

n = V2^^£ = -e-^°P = -e-^°- . (4.3) 

As we discussed in the previous section, when we integrate the equations of motion nu¬ 
merically we fix X = 1 at the horizon, and use the rescaling (3.15a) to meet the boundary 
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Figure 2; (Color online) Values of the condensate as a function of the temperature for = 

0 and qL = 5/2. Each curve corresponds to a different value of (see inset). Dashed 

(solid) lines indicate first (second) order phase transitions. 


condition x(r —oo) = 0. Applying this rescaling, the temperature and the chemical potential 
of the dual theory are given by 

T = , /i = , (4.4) 

in terms of T, p,, and xo defined in eqs. (3.17), (4.2c) and (4.2b) respectively. We then 
use the chemical potential to construct dimensionless ratios of the physical quantities in our 
problem, and present our results in terms of these. For example, the corresponding ratios for 
the temperature and free energy read 

T _ T 2k2 D _ —2x0 ^ 

/r p ' V2 ^ p'^ ' 

where the left-hand side of these expressions corresponds to field theory quantities and the 
right-hand side to parameters from the numeric solution. 

From equation (3.22) we obtain that the condensate density of the marginal operator can 
be read from the UV asymptotic solution as 

2K^^ = 3nT3. (4.6) 

In figure 2 we plot the value of the condensate (actually rs, which is related to the con¬ 
densate by means of eq. (4.6)) as a function of the temperature for a single value of the scalar 
charge qL = 5/2, and different values of the tension, Tb. The outer lines correspond to phase 
transitions with larger critical temperature and vice versa. For the solid lines a fixed value 
of the temperature corresponds unambiguously to a value of the scalar condensate. At the 
critical temperature the value of the condensate is turned on continuously, and the phase 
transition is of second order. For the dashed lines there is a region around the critical tem¬ 
perature in which metastable configurations exist for a given temperatnre. For T > Tc the 
thermodynamically preferred phase is the one with no condensate, and for T < Tc the one 


(4.5) 
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Figure 3: Plots of the the free energy in the normal (black, dashed line) and condensed phase 
(red, solid line) as a function of the temperature for m? = 0 and q L = 5/2. 
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Figure 4: Phase diagram for m?‘L?‘ = 0 and qL = 5/2. Blue disks correspond to second order 
phase transitions and red boxes to first order ones. The inset is a zoom on the rightmost 
points. Beyond = 16/3 there is no condensed phase at any temperature. 


with the largest condensate. There is a finite jump in the value of (O 3 ) as the temperature is 
lowered, and the phase transition is first order. 

The value of the critical temperature is determined by studying the free energy of the 
system, and in figure 3 we present two examples of this calculation. In both graphs the line 
corresponding to the free energy of the normal phase was obtained from eq.(3.2) as = —P. 
The left-hand side graph, figure 3a, corresponds to a second order transition, while the plot 
on the right-hand side, figure 3b, corresponds to a first order one: the derivative of the free 
energy with respect to the temperature presents a finite jump. 

In figure 4 we plot the value of the critical temperature as a function of the tension Tb 
for the same value of the charge as before, namely qL = b/2. To ease the identification 
of the phase transition order for a given tension, we have indicated the second (first) order 
phase transitions with blue circles (red squares). As is clear from the plot, for this value 
of the charge there exists a region, corresponding to intermediate values of the tension, for 
which the phase transition is of the first order kind. However, this region does not extend to 
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qL 

Figure 5: Order of the phase transition in the qL - ThL?' plane. The red, shadowed region 
denotes values of the parameters where a first order phase transition occurs, and the blue 
striped region surrounding it to values where a second order phase transition does. The thick 
solid and dashed lines are the curves in figure lb. In the white region there is no condensed 
phase. 


arbitrary values of the charge qL, as one can see in figure 5. In this plot we present a scan 
of the parameter space {Th,q), showing that only for a closed region in that plane do first 
order phase transitions occur. At zero values of the tension, or for values of the parameters 
where the critical temperature is infinitesimally small, the phase transition is always second 
order. Pictorially, the red, shadowed region in hgure 5 is fully surrounded by values of the 
parameters where the transition is second order. 

An interesting observation about the plot in figure 5 is that the solid thick line separat¬ 
ing the region with no condensate (large values of the tension) and the region with second 
order phase transitions (striped region) corresponds to an holographic BKT quantum phase 
transition governed by the instability occurring at the AdS 2 throat described in section 3.2. 
Indeed, for the marginal operator under consideration this instability takes place along the 
curve 




(4.7) 


where we have defined L^ as the tension that saturates the AdS 2 BF bound at extremality 
as a function of the charge of the scalar field. On the other hand, the dashed line in figure 
5 corresponds to the instability beyond the AdS 2 throat, also discussed in section 3.2, and 
describes a second order quantum phase transition. These two lines (solid and dashed) are 
the same ones as those in figure 1. To understand better these phase transitions it is worth 
studying in detail the behavior of the system at zero temperature. Although this study will 
be the subject of the next section, let us point out that in the inset of figure 4 it is shown 
how the critical temperature approaches zero exponentially fast when —)• = 16/3, 

following the behavior in equation (3.II). 

We finish this section with some representative plots of the behavior of the entropy density 
as a function of the temperature in figure 6. Both in the normal and broken phases the entropy 
density is given by the horizon radius squared (r^) up to some proportionality constants (see 
eq. (3.2b)). In our numerical integration we keep = I by making use of the scaling 
symmetry (3.15b), hence the rescaling invariant ratio s//i^ is given by 


s K 
2ti 


o-xo 


(4.8) 
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Figure 6: Plots of entropy density (in units of k?/{2t^)) as a function of the temperature at 
m? L?' =0, qL = 2.5 and two different values of the tension corresponding to (a) 1®* order 
and (b) 2”'^ order phase transitions. 


where fi and xo respectively the asymptotic values of At and x defined in eq. (4.2). In 
the plots in figure 6 the dashed, black line corresponds to the normal phase, and it goes to a 
constant in the zero temperature limit. The continuous, red line corresponds to the condensed 
phase; near T = 0 it vanishes quadratically with the temperature. 


4.1 Quantum phase transitions 


In this section we will study the zero temperature limit of the system with a massless scalar, 
which will allow us to characterize better the quantum phase transitions described above. 

The zero temperature geometry corresponding to the ground state of the condensed phase 
consists of a domain wall interpolating between two AdS^ with the same radius. We construct 
it by integrating numerically from the IR (r —)■ 0) towards the UV(r —)■ oo). In the IR the 
metric is aAdS^ with the same AdS radius as in the UV, constant scalar field, and vanishing 
gange field. Conformal symmetry is therefore recovered in the IR end point of the RG flow for 
the condensed phase.® The asymptotic solution around the origin of spacetime is, for q ^ 0 
(recall that we are setting L = 1), 


S' ~ ^1 — (a + I) + O ^ 


X-Xo 


q'XO + 1)(Q? + 2) ^2a _|_ Q 

4 a VI’ 


(r2“ + 0(r^“)) 


'j- 


■\/{a + l)(a + 2) q^ e^° 


V2q 


4 a (3 + 2a) 


r2“ + 0(r^")^ , 


(4.9a) 

(4.9b) 

(4.9c) 

(4.9d) 


where a > 0 is a free parameter. 

In (4.9) we have imposed four conditions by choosing the IR behavior of the fields, and 
we need to set two further boundary conditions in the UV (r —)• oo). These are as before 
X = 0, and the vanishing of the source term in the UV asymptotics of the scalar. Again, we 
can make use of the scaling symmetry (3.15a) to make x 0 as r —)• oo. Then, by tuning 

®Tlie same phenomenon occurs in the Einstein-Maxwell-scalar model of [28], whose treatment we follow 
closely in this section (see also [33]). 
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Figure 7: Plots of the functions in our ansatz as a function of r j at zero temperature, 
re? I? = 0, q L = 2. The different curves correspond to tensions = 1/2 (red, dotted), 

= 3/2 (cyan, dashed) and = 7/3 (blue, solid). 


the constant a in our IR solutions we ensure the vanishing of the source term in the UV of 
the scalar. 

In the probe limit Tj, = 0, the equations for the metric can be integrated readily as 
g = and x = Xo- Then, numeric solutions for At and r, that depend on the value of the 
charge of the scalar, correspond to a probe spacetime filling brane in an AdS^ spacetime, 
with action given by just the square root term in (2.7). For finite values of the tension we 
integrate numerically the whole set of equations (2.15) using the IR asymptotic solution (4.9) 
(expanded up to a higher order in r^“) as a seed. 

In figure 7 three examples of the numerical integration are shown for a charge qL = 2 
(in each case a has been tuned to a value such that the source term of the scalar is |ro(r = 
10^)1 < 10“^^). In view of these results we shall next study how the system behaves as one 
approaches the critical values of the parameters for which the instabilities studied in section 
3.2 occur. 

First, remember that for vanishing tension the system becomes a probe in AdS 4 ,^ and 
indeed at low values of Tt, the value of g/r'^ is close to the AdS^ value gjr^ ~ 1 everywhere. 
As the tension is increased a dip appears in g{r\ with the minimum of the dip having a positive 
value. We know from the analysis in section 3.2 that there exists a critical value of the tension 

beyond which the condensed phase does not exist anymore. When this critical tension is 
reached, Tf, L^, the minimum of g occurs at [g/r'^)\r=rQ = 0, signaling the appearance 

of an extremal black hole with horizon radius tq. Moreover, this numerical solution coincides 
accurately with the extremal black hole normal phase solution, given in eqs. (3.1), for r > rg. 

Regarding the scalar r, for the critical case Tj,L^ = L? the scalar is different from 
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{h) qL = 3/2 


Figure 8: Condensate at zero temperature for = 0. The inset in the left-hand side plot 

is a zoom on the exponential tale. The dashed line corresponds to expression (3.11) for the 
holographic BKT phase transition with a fitted value of Ai/j. The solid line connects linearly 
the numerical data. 


zero only inside the horizon, r < tq, whereas for lower values of the tension its profile is non 
trivial everywhere, indicating (a) the existence of a condensate for Tb related to 

the normalizable mode of the scalar at the boundary, and (b) that in the condensed phase the 
charge at zero temperature is extended in the bulk of the spacetime, whereas for the normal 
phase it is contained in the extremal black hole [33, 28]. 

The metric function x behaves similarly to the scalar field r (notice that in figure 7 we have 
scaled its value at the origin to one for presentational purposes), being zero outside the horizon 
for the critical case, and extending into the bulk for lower values of the tension. Generically, 
the value x(0) ^.t the origin grows with the tension up to a maximum at (and 

it is zero for = 0). 

Finally, the gauge field profile vanishes for r < rq when TbLp' = L?', and extends down 

to r = 0, vanishing there, for lower values of the tension. 

As a further check of our numerics, we have verified that the values for the condensate 
obtained at zero temperature coincide with the values one gets by extrapolating the low tem¬ 
perature behavior observed, for example, in figure 2. Actually, the whole numerical solution 
at finite temperature asymptotes to the T = 0 one as one reduces the temperature. 

To conclude this section let us corroborate that, as stated in section 3.2, when the quantum 
phase transition occurs due to the violation of the BF bound of the near horizon AdS 2 
spacetime, this phase transition is a holographic BKT (continuous) transition. While this is 
not the case for the second kind of instability discussed in section 3.2 

In a BKT phase transition, the behavior of the condensate at zero temperature as the 
critical point TbL? = is approached must be given by the expression (3.11). We shall 

then compute the value of the condensate at T = 0 as a function of the tension, in the 
vicinity of the critical point (4.7). We present our results in figure 8a, where on the inset 
we provide a zoomed-in version of the exponential tale. The dashed line there corresponds 
to the prediction of eq. (3.11) after a one parameter fit using the position of the rightmost 
point shown in the plot. Instead, in figure 8b we provide the same plot for a value of the 
charge for which the instability at zero temperature occurs away from the horizon (he., it 
is the second kind of instability discussed in section 3.2). It is clear from this figure that 
there is no exponential holographic BKT behavior. This was expected since the instability is 
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not caused by the violation of the BF bound in AdS 2 , and the quantum phase transition is 
actually second order. 


5 Phase diagram for the model with a A = 2 operator 


In this section in order to study the theory where the condensing operator has dimension 
A = 2 we take the mass of the scalar field to be m? = —2. With this choice the scalar field 
behaves asymptotically as 

r~A_)_A_| -^ (|7-_).oo) (5.1) 

Both Ti and T2 are normalizable modes, and hence two different quantizations are possible [34], 
The two theories correspond to the cases where ri {T 2 ) is dual to the source of a dimension 
A = 2 (A = 1) operator and T 2 (ri) is dual to the VEV. 

For the remaining functions in our ansatz, the equations of motion imply that in the UV 
(r —00) 


g r 


. Tb rf £ 

1 H-H-h • 

2 7-2 7.3 ' 

Th ATt, T1T2 


X ^ Xo + —^ + 


A, 


_2 

- Q 

g -h 

r 


+ 


(5.2a) 

(5.2b) 

(5.2c) 


From now on we pick the quantization in which the operator has dimension A = 2, and 
correspondingly ti is interpreted as the source. For the scalar to condense spontaneously we 
impose ri = 0. In this case the expressions for the physical quantities given in eqs. (4.3) and 
(4.5) hold, with the value of £ to be read from the UV asymptotics of g given in eq. (5.2a). 
The condensate density of the dual operator obtained from (3.22) reads 


2K^^ = nT2. (5.3) 

The same procedure as described in the previous section is employed to construct the 
numerical solutions corresponding to the condensed phase of the system at finite temperature. 
In order to illustrate these solutions, a plot of the condensate versus the temperature is shown 
in figure 9, while for the phase diagram in the q L ~ T},LP‘ parameter space we refer the 
reader to figure la. Recall that in that figure the solid line corresponds to holographic BKT 
transitions at zero temperature, whereas the dashed one corresponds to T = 0 second order 
phase transitions. Moreover, we have found that, as opposed to the case in the previous 
section, for a relevant operator of dimension A = 2 the phase transition to the condensed 
phase is always of second order at finite temperature. 

For the case at hand, ie. m? = —2, we have not constructed the solutions correspond¬ 
ing to the ground state (zero temperature) of the condensed phase. However, by studying the 
behavior of our solutions at very low temperature one can gain intuition on the nature of the 
ground state. The result is presented in figure 10. First, notice from the plot of the function 
gjr^ that the radius of the horizon (where giji^ = 0) decreases, in units of the chemical po¬ 
tential, as the temperature is decreased. Precisely at the horizon At{rh) = 0, while and 

r(r/i) go to constants; as required by the boundary conditions. However, the most interesting 
feature of these plots is the region that appears at small radius in the T —^ 0 limit. In this 
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Figure 9: (Color online) Values of the condensate as a function of the temperature for 
rm? = —2 and qL = 2. Each curve corresponds to a different value of the tension (see 
inset). We have checked that all the phase transitions are of second order. 


region both g/r'^ and x stabilize to a constant, hinting to an AdS^ geometry for the ground 
state. The scalar field r diverges towards the horizon, while the gauge field At is suppressed 
exponentially as r —)■ 0. This picture is confirmed by the equations of motion (2.17) in the 
small radius limit. Indeed, at leading order in the IR, one can find a solution where the 
matter fields r and At give exponentially suppressed corrections to the geometry. The metric 
becomes that of AdS^^ with radius (1 + Th L'^/Q)L, while the scalar diverges as 7--4/(6+r6L ) 

when r —)• 0, whereas the gauge field vanishes exponentially At ~ exp r . yVe 

leave for the future the construction of the zero temperature solution for the whole range of 
the radial coordinate, which is expected to be a domain wall interpolating between two AdS^ 
geometries [35]. Finally, notice that when r//r —)• oo, the gradient of r approaches a constant 
that remains basically unchanged for many (low) values of the temperature, corresponding to 
the flat parts of the curves in figure 9. 

Regarding the entropy density, at low temperatures this quantity behaves as dictated by 
the AdSn geometry, namely it vanishes quadratically with the temperature. 

To finish this section let us comment on the region for which the condensed phase does not 
exist for intermediate values of the tension Tf, (as it can be seen in figure la this happens 
for values of the charge above qL = lj2 and below qL ^ 1.936). From the discussion above, 
the zero temperature transition that occurs at small values of the tension (along the dashed 
line of Fig. la) is of second order, whereas the T = 0 transition at large values of the tension 
(continuous line in Fig. la) is of the holographic BKT type. Moreover, the typical critical 
temperature for large values of the tension is several orders of magnitude smaller than the 
chemical potential, and it grows slowly as the tension is increased (this last fact is true also 
for values of the charge such that the condensed phase exists for any value of the tension, 
namely those above qL ^ 1.936). 
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and Tfo = 1. The different curves correspond to decreasing temperatures from T = Tc to 
T = 0.015Tc (thicker, orange line). 


6 DC and optical conductivities 


From the phenomenological point of view it is of great interest to calculate the retarded 
current-current correlators, and in particular the conductivities associated to the different 
currents. Holographically we achieve this by studying the set of fluctuation modes that couple 
to the spatial perturbation of U{1) gauge fields. In the UV we require our conhguration to 
realize an electric field of constant modulus and frequency uj along the x direction, while 
in the IR we impose ingoing boundary conditions at the horizon [36]. The conductivity is 
then computed as the quotient between the normalizable and non-normalizable modes of the 
spatial, zero momentum perturbation of the U{1) field, divided by iu. 

Recall that, as discussed in section 2, in our setup there are two gauge fields corresponding 
to the axial and vectorial U{l)s. In the absence of chemical potential along the U{l)v the 
fluctuations corresponding to switching on an electric field in either of the U{l)s do not mix 
with each other, and hence we can study them separately. We start the discussion with the 
vectorial conductivity. Therefore we consider the following perturbation of the U{l)v gauge 
field: 


6V = e ^‘^^Vx{r)dx . 


( 6 . 1 ) 


To write the Lagrangian quadratic in the fluctuations we must consider the most general 
version of the model, for which the DBI term is given by eq. (2.1). The fluctuation (6.1) 
decouples from any other mode, and the relevant piece of the action is then given by 


/’{ 2 ) 


r\j 


Tb 1 


V-det , 


( 6 . 2 ) 
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with A4 and W defined in eqs. (2.13 -2.14), and the field strength associated to the [/ (l)y 
field. The effective, radial-dependent coupling reads 


1 


V(t) 




(6.3) 


and is proportional to the potential and the quotient of the metric and open string metric 
determinants. The resulting equation of motion reads 


v" -I- dr log 


lats 


^ y^rr 


Vx. — 0 . 


(6.4) 


Since translation invariance is unbroken in our model we expect that charge will not dissipate 
in time, leading to a delta peak at zero frequency in the conductivity (equivalently, a l/w pole 
in the imaginary part). However, in the Lagrangian (6.2) there is no mass term, which implies 
the absence of such a delta peak. The reason for this is that the charge density associated 
to the U{l)v field is exactly zero in our setup,and therefore the electric current, which is 
then due to pair creation only, carries no net momentum. 

We will as usual read the AC conductivity from the asymptotic behavior of Vx{r). For 
large r the solution of (6.4) behaves as 


( 1 ) 


Vx = H-h o(r ^) . 


(6.5) 


We will solve (6.4) numerically imposing ingoing boundary conditions at the horizon, and 
then read the conductivity as 

vW 

^-(oT ■ 

VjJ UW; 

Notice that the numerically obtained quantity a coincides with the conductivity up to a 
normalization 


a 


V, 


LO = 


A 

2^2 


o- w 


(6.7) 


In particular lim(^_).oo a = 1. 

As for the DC conductivity in the U (l)y sector, isotropy and the effective Maxwell action 
(6.2) allow us to follow the procedure of [37] and express the DC conductivity in terms of 
background functions evaluated at the horizon. Both for the normal and condensed phases 
the result can be written as 


y n V{r+) 


( 6 . 8 ) 


where r+ is the value of the scalar at the horizon r = r_|_. Notice that this expression becomes 
even simpler in the normal phase where the scalar vanishes and therefore H(t+) = 1 

Let us see what we can learn about in the condensed phase for the two cases analyzed 
in this paper. In the model with m? L?' = 0, i.e., that with a marginal operator, the potential 
is just F = 1, and thus is determined solely by the value of A[ at the horizon. Notice 
that A((r+) vanishes in the zero temperature limit (see the IR expansion (4.9c)), and hence 
at T = 0 the DC conductivity is a non zero constant. 

For the model with L? = —2, the DC conductivity of the U{l)v current takes the form 


— 

^DC — 


.-2Ti 


2K^^l-A'{r+y ■ 


(6.9) 


^°Notice that in all our results in this section the chemical potential will be always that along the U{1)a- 
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Figure 11; DC conductivity of the vectorial current for 


-2,qL = 2, andTfoL^ = 0.25. 


Interestingly, as shown in figure 10, as T —)• 0, t+ diverges while A^{r^) —)• 0, hence aY)c is 
suppressed as . Therefore, in this sector (for which the charge density is zero), the system 
behaves as an unconventional insulator.We provide an example of this feature in figure 11, 
where we observe that below a certain temperature the DC conductivity is suppressed. This 
suppression is an effect of the effective metric that governs the fluctuations of the vectorial 
field (see eq. (6.2)). This effective metric is sensitive to the value of the condensed scalar, and 
even when we have not broken the U{l)v symmetry we observe a pseudogap (the conductivity 
is small but not exactly zero) in the corresponding conductivity. This is a direct effect of our 
DBI model. 

We shall now focus on the conductivity associated to the axial U{1)a held. We must then 
consider the perturbation 

6A = e~^‘^^ax{r)dx , (6.10) 

which is, as usual, coupled to the fluctuation 6gtx of the metric due to the finite chemical 
potential associated to the U{1)a symmetry. However, a constraint relating drdgtx to ax, 
arising from the Einstein equations, allows us to decouple the ax mode completely and solve 
a simple linear differential equation of second order.The effective Lagrangian for the axial 
perturbation is then given by (6.2) with the substitution Fy —)• Fa and the addition of a mass 
term. The resulting equation of motion for ax{r) takes the form 


a" + dr log 


1 ,_ 1 yatt 

-^V-detMW^^W a'x-uj^^ax + 


Y{r) 


" ' 4 (e-x 5 - 2g2 t2) 

with the mass term given by the following expression 


fla; = 0 , (6.11) 


Y = 2q^ 


re ^ {l + 2gT'^ - {Af + 2q^ A^ g ^ + 2gT'^)))+2 t A'^A^^dr 


+ nv{T) 
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/2 


^1 + 2 5 r'2 — [A'^ + 2q‘^ A^ g~^ {l + 2g r'^)) 


Ty/AlJ 

( 6 . 12 ) 


^^Note that a similar feature was displayed in [38] where in the absence of charge density the model exhibited 
an exactly vanishing DC conductivity (see also [39]). 

^^This is not unfamiliar in the contest of holographic superfluids, see for example [5] for explicit expressions 
in a simpler setup. 
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Figure 12: Real part of the (a) vectorial conductivity and (b) axial conductivity in the 
condensed phase as a function of the frequency for m? Lp' = 0, q L = 2.5 and = 0.15. 

The line in the plot is darker the lower the temperature is, taking values from T/T^ = 
{1, 0.81, 0.41, 0.12} (the two last values are almost indistinguishable in the plot). 


The UV expansion of ax is of the same form as (6.5) and thus to compute the conductivity 
we solve the equation (6.11) numerically with infalling boundary conditions at the horizon, 
and read the conductivity using eq. (6.6). Notice that Y(r) now implies the presence of a 
delta peak at zero frequency in the U{1)a conductivity. Whenever we talk about the DC 
conductivity in this case we refer to the finite part that can be read from the zero frequency 
limit of the optical conductivity. 

Let us now present and discuss the plots with our numerical results for the optical conduc¬ 
tivities. First, in figure 12b we plot the real part of the axial conductivity for the condensed 
phase of the model with m = 0. These results are very similar to what was obtained for the 
conductivity of the minimal holographic superconductor in [5]. Namely, the AC conductivity 
approaches a constant at large frequencies (as determined by the AdS^ asymptotics of the 
solution), while at low frequency a pseudogap appears when the temperature is low enough. 
In this hgure we also observe that the AC axial conductivities for temperatures smaller than 
T ~ O.dTc remain basically invariant. The reason is that for low values of the tension, and 
temperatures below ~ 0.4Tc, the dynamics is well approximated by that of the zero temper¬ 
ature solution discussed in section 4.1 (whose geometry corresponds to an AdS domain wall) 
This can be argued by looking at figure 2 where one can see that at T ~ OATc the condensate 
has already stabilized at the corresponding zero temperature value. In figure 12a we plot 
the real part of the vectorial conductivity also for the condensed phase of the model with 
m = 0. As discussed above there is no delta peak in this sector, and moreover no pseudogap 
appears either. Again the conductivity practically reaches the zero temperature curve for 
temperatures below ~ OATc. Moreover as T —)• 0, the DC conductivity tends to the AdS 4 
conformal value Tb/{2K^), as predicted by eq. (6.8). Notice that in figure 12a the two lines 
with larger temperature go to a constant at zero frequency, although this is not visible since 
we have truncated the axis for presentational purposes. 

In figure 13a we provide the U{l)v optical conductivity in the condensed phase for the 
case with mP Lp = —2. Let us first emphasize again that no delta peak appears in the 
zero frequency limit. Indeed, as we discussed below eq. (6.9) the DC conductivity is finite, 
and moreover is highly suppressed as the temperature approaches zero. Next, also for low 
temperature, a pseudogap is clearly visible at low frequencies. The existence of this pseudogap 
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Figure 13: Real part of the (a) vectorial conductivity and (b) axial conductivity in the 
condensed phase as a function of the frequency for m? L'^ = —2, qL = 2 and T^L'^ = 0.25. 
The line in the plot is darker the lower the temperature is, taking values from T/Tc = 1 
(lighter) to T/Tc = 0.09 (darker). 


is a consequence of the behavior of the tachyon in the IR. As we have seen in figure 10, the 
tachyon diverges towards the horizon, and hence the factor V{t) appearing in the effective 
coupling (6.3) has the effect of a soft wall suppressing the spectral density at low frequencies 
as in [40]. Since the AC conductivity satisfies a sum rule, the suppression of the a; = 0 
conductivity has to be compensated by an increase of the spectral density at finite values of 
the frequency. This shift of spectral weight is confirmed by the numerical solution; at low 
temperatures a peak forms at the end of the pseudogap. After that peak there is a region with 
approximately constant conductivity, followed by a decrease towards the value determined by 
the asymptotic AdS^ geometry at large frequencies. 

In figure 13b we present the real part of the AC conductivity for the axial current. As 
discussed above, a delta peak at zero frequency is present both for the condensed and normal 
phases. This delta peak is visible via the dispersion relations through the associated 1/uj 
pole in the imaginary part of the conductivity. One can conhrm that the weight of the delta 
function receives a new contribution (proportional to the superfluid density) once the system 
is in the condensed phase. Indeed, as shown in figure 14 (where we study the evolution of the 
residue of the pole of the imaginary part of the conductivity as a function of temperature) one 
observes a continuous change in the behavior of the zero-frequency pole as the system enters 
the superfluid phase. We also observe that at low frequencies and low enough temperature 
the AC conductivity presents a pseudogap. Notice that the width of this region is largely 
independent of the temperature, Ug ~ 0.42/i. 

To end this section, let us insist on a feature of hgure 13. At very low temperatures (about 
10% the critical temperature) a sharp peak becomes visible at the end of the pseudogap 
region. This peak becomes higher and narrower as the temperature is decreased, and we have 
checked that it appears quite generically when varying parameters of the theory, being easier 
to observe for low values of the tension Tf, Lp'. We comment on its possible significance in the 
next section. 
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Figure 14: Pole of the imaginary part of the axial conductivity Ug for the m? = —2 case. 
Us results from a fit of Im[(T] to the function ns/{io/n), for values of a;//i < 0.1. The solid red 
line corresponds to the normal phase, while the dashed black line to the broken phase. 


7 Discussion 

In this work we have constructed a holographic superfluid with matter transforming in the 
fundamental representation of the gauge group. The main novelty of our model comes from 
the use of the tachyonic action describing a spacetime filling D3-D3 in an asymptotically 
AdSi BH geometry. The setup contains a scalar field, the tachyon, charged under the U (1)a 
included in the global U{l)xU{l) symmetry supported by the D3-D3, but neutral under the 
diagonal U{l)v- We have turned on a finite chemical potential corresponding to the U{1)a, 
and found that below some critical temperature the tachyon condenses breaking U{l)xU{l) —>■ 
17(l)y, and hence realizing a superfluid phase transition. We have focused on two classes of 
theories, one where the operator dual to the tachyon is marginal (A = 0), and the other one 
corresponding to a relevant operator with A = 2. 

We have shown that there are two gravitational types of instability that drive the tran¬ 
sition to a condensed phase, one of them giving rise to a holographic BKT phase transition 
at zero temperature, whereas the other leads to a second order one. For the case with a 
marginal operator we have constructed the zero temperature solution and shown that the low 
temperature dynamics is governed in the IR by a conformal fixed point with the same central 
charge as the fixed point governing the UV. 

Once the mass of the tachyon is fixed, there are two further parameters in the model; these 
are the tension of the branes (which controls the amount of backreaction), and the charge of 
the tachyon. We have studied the phase diagram of the system as a function of these two 
parameters plus the temperature. When the operator that condenses is marginal there are 
two different types of finite temperature phase transitions. Depending on the explicit values 
of the charge and tension, the finite temperature transition can be first or second order. In 
the former case there exist metastable phases close to the critical point, whereas in the latter 
the value of the condensate provides an order parameter with mean field exponent. When 
the operator that condenses has dimension A = 2 the phase transition at finite temperature 
is always second order. 

In the last part of this work we have studied the AC and DC conductivities for our setup, 
both for the current associated to the U{l)v and the U{1)a symmetries. In the U{1)a, the 
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conductivity shows features very similar to those of the minimal holographic superconductor 
of ref. [5]. It displays a delta peak at zero frequency both in the normal and condensed phases, 
with the weight of the delta function receiving a new contribution when the system enters 
the broken phase. As expected, the axial optical conductivity presents a pseudogap at low 
temperatures. The U{l)v DC conductivity is instead finite both in the normal and broken 
phases. This is due to the absence of charge in that sector, for which only the pair produced 
charged carriers contribute to the conductivity. We have produced a neat expression for this 
DC conductivity in terms of the horizon data, showing that, for the case where the operator 
dual to the tachyon has A = 2, as the temperature is lowered the DC conductivity goes to zero 
faster than exponentially, and thus the system behaves similarly to an insulator in this sector. 
As for the optical conductivity of the U{l)v sector, a pseudogap appears at low temperatures 
for the theory with a A = 2 operator due to the scalar potential acting as a soft wall in the 
IR for the fluctuations. As expected, no pseudogap is observed for the theory with a marginal 
operator (where the tachyon is massless, and therefore the potential is trivial). It would be 
interesting to investigate what are the effects of a non trivial U{l)y chemical potential on the 
conductivities. 

Finally, let us comment on an interesting feature of the optical conductivities for the model 
with a A = 2 operator. At low temperature, and for low values of the tension, there appears 
a sharp peak at the end of the pseudogap region. This peak is the signature of a quasinormal 
mode (QNM, a normalizable solution of the fluctuation equation) located at a specific position 
in the complex frequency plane, D. For larger temperatures the position of this QNM has 
a negative imaginary part Im [D] ~ —T. When the temperature is lowered this mode comes 
closer to the real frequency axis, and a peak with a certain width and height appears in the 
spectral function; a quasiparticle. If eventually the mode becomes real, Im [If] = 0, the peak 
becomes a delta function located at cj = D, and corresponds to a particle of the spectrum 
of the theory with mass M = D. It is then natural to conjecture that the theory with a 
relevant operator of dimension A = 2 has, at zero temperature in the condensed phase, a 
particle precisely at M = fl. In principle more massive modes would be present as well, and 
be observable in the spectral function at extremely low temperatures, but we were not able to 
obtain enough control on the numerics to observe them. We plan to come back to this point 
in future work [35]. 

A second line of future research would consist in the study of unbalanced superconductors. 
In order to do so one needs switch on a chemical potential along the diagonal U{l)v- Since 
the tachyon is not charged under this U{\)v one can interpret this second chemical potential 
as an imbalance of populations in a mixture of two species [17]-[19]. This is an interesting 
line of study, since at weak coupling the so called LOFF inhomogeneous condensed phase [41] 
is expected to appear [42]. In this context, the appearance of first and second order phase 
transitions in the phase diagram of the theory with a marginal operator is already hopeful. 
As reviewed in [18], in the phase diagram of unbalanced superconductors there is a region, 
at large values of the imbalance, where the phase transition becomes first order, and it is in 
this region where the LOFF phase is predicted to appear. It would then be interesting to 
consider the model with A = 0 and values of the parameters close to the region in the phase 
diagram where the phase transition becomes first order (see figure 5), and study the effect of 
a non zero chemical potential along the diagonal U{\)v- Notice that while the two D(l)s in 
[18] did only interact with each other through their backreaction on the geometry, thanks to 
the non-linearity of the DBI action, interacting terms will be present in our setup. 

Finally, a more ambitious continuation of this work would be to consider generalizations 
of this model along the lines of [43] with a non trivial dilaton and consequently more general 
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potentials depending both on the tachyon and the dilaton. This approach could lead to an 
interesting phenomenology as in [44] for Einstein-Maxwell theories. 
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